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I. INTRODUCTION 

An important way to understand the interactions 
among particles is via the study of scattering. In related 
experiments, the initial status of a scattering state is re- 
quired to be precisely controlled. This is relatively easy 
for charged incident particles because their initial mo- 
mentum can be tuned by adjusting the electromagnetic 
forces imposing on them. For neutral incident particles, 
the initial momentum is in general difficult to control pre- 
cisely. However, the recent progress in the technology of 
trapping neutral atoms via optical trap might open a new 
way for studying the scattering of neutral particles with 
precisely controllable initial status [IH1]- In this paper 
a model of a device is proposed and related theoretical 
calculation is performed to show how the mentioned scat- 
tering is realized. It turns out that in this device, as we 
shall see, the collisions among particles occur regularly 
and repeatedly. Thereby the effect of each individual 
collision can be accumulated. This would be helpful for 
the understanding of the very weak interactions among 
neutral atoms (molecules). 

Traditionally, the scatterings would at most have two 
incident channels (say, in the experiments with head-on 
colliders). However, in the following device, three or 
more incident channels can be realized. As an exam- 
ple, a three-body scattering with three incident channels 
is chosen to be studied. This is a generalization of a pre- 
vious work on a two-body scattering in a trap with two 
incident channels Q. 

It is assumed that, in the beginning, there are three 
narrow optical traps located at the three vertexes of a 
regular triangle, and each optical trap provides a har- 
monic confinement. The total potential is U p (r) = 

TjMujp Ylj=i l( r — a j)| 2 ' where &j points from the origin 
to the j — th vertexes, and M is the mass of a parti- 
cle. When the center of the triangle is placed at the 
origin, = a. It is further assumed that each trap 
contains an atom in the ground state of a harmonic os- 
cillator, the three atoms are identical bosons with spin 
zero, and uj p is large enough so that the atoms are well 
localized initially. Suddenly the three narrow traps are 
cancelled. Instead, a broader new trap located at the 



^Muj 2 r 2 is created, ui < uj p . Then, the 
The evolution is affected not 



origin U evo i(r) 
system begins to evolve 
only by U evo i(r) but also by the atom-atom interaction 
V(|rj— Tj |). In what follows the details of the evolution is 
studied, three-body head-on collisions occurring repeat- 
edly are found, and the effect of interaction is demon- 
strated. 



II. INITIAL STATE 



We shall use Hui and y/h/Moo as units of energy and 
length. The symmetrized and normalized initial state 
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where implies a summation over the permutations 
P1P2P3, and 77 — uj p /uj. Without loss of generality, a 3 
is given lying along the Z-axis, while the triangle is 
given lying on the X-Z plane. For convenience, three 
sets of Jacobi coordinates denoted by a, /3, and 7, re- 
spectively, are introduced. The coordinates of the a set 
are defined as, r = r 2 — ri, R = r 3 — (ri + r 2 )/2, and 
R c = (ri+i"2+r3)/3. The other two sets can be obtained 
from the a set by cyclic permutations. In terms of the 
a set, we introduce the harmonic oscillator (h.o.) states 

0niln( s ) = fnl( s ) Y lm(s) as basis functions, where s = r, 
R, or R, 



tonian — ^7 V| 



They are normalized eigenstates of the Hamil- 
\nS 2 with the eigenenergy 2n + 1 + 3/2 



and with the angular momentum / and its Z-component 
m, where fi = 1/2, 2/3, and 3 when s = r, R, and R c , 
respectively. Then the initial state can be expanded as 
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where the first factor is for the cm. motion which is 
completely separated from the internal motion, 
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The notation q denotes a set of quantum numbers n, I, N 
and L, and II = (— l) l+L is the parity, 



$ Jrn n,(r,R) = [^2 /2) (r)^(R) Jm ] 
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where I and L are coupled to J and m. 

GjmUq — 'y^^a n 'l'mbN'L'Cl>m,L'o(l + ) 

x [V + Af 7 (/3 -> a) + ^'' 7 ( 7 -> a)] (5) 
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Q 7 ™ L , is the Clebsch-Gordan coefficients, q' is for the 
set n', I', N' and L' to be summed up, 8 r and 4> r are the 
spherical polar coordinates of r, and so on. 
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is the bracket of transformation between the j3— and 
a— sets (the superscript a is usually ignored), which is 
called the Talmi-Moshinsky (T-M) coefficients. Their an- 
alytical expression can be found in 0-0] ■ In general, the 
symmetrization would cause the appearance of all three 
set of coordinates. However, by using the T-M coeffi- 
cients, Eq. ([2]) contains only the a— set so as to facil- 
itate greatly the calculation. Incidentally, since ^1 is 
symmetrized, / (included in q) of Eq. @ must be even 
and therefore (— 1) L = II. In principle, the right side of 
Eq. (|5J) should contain infinite terms. However, the over- 
lap between VP/ and higher h.o. states are very small. 
Say, if, n, N, N c are all smaller than 12 and I, L, L c are 
all smaller than 23, with the parameters specified below, 
the overlap of the right side of Eq. (|5J) with itself is equal 
to 0.99997. This implies that higher h.o. states can be 
safely ignored. 

III. HAMILTONIAN AND ITS EIGENSTATES 

The evolution is governed by the Hamiltonian H evo i 
containing U evo i (r) and the interaction. When the Jacobi 



coordinates are used it can be separated as 

H evo i — H c + H in (9) 

where H c = — ^Vp^ + ^R 2 describes the cm. motion, 
and 

Hin = -V? + \r 2 - \vl + \ R2 + T, V(\*i ~ (10) 

describes the internal motion. Since the eigenstates of H c 
are well known, if the eigenstates of Hi n are also know, 
the evolution starting from any initial state can be un- 
derstood. In what follows, the symmetrized eigenstates 
are obtained via a diagonalization of Hi n in a limited 
space. In this way, only approximate solutions can be 
obtained. Then, we increase the dimension of the space 
until a better convergency is achieved. 

When $,/ m n 9 (r,R) as defined in Eq. (gj) are used as 
basis functions (where I is restricted to be even as men- 
tioned), the matrix elements of Hi n is 



( $j'm'nv(r,R)| H in | $j mrig (r,R) ) 

= Sj>jS m : m Su'll [ 5q'q(2n + 2N + I + L + 3) + Sl'l^N'N^L'LVn'nl 
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where q" denotes the set (n"l"N"L"), and the implica- tion of q'" is alike, 

Vwm = J r 2 dr f%' 2 \r) V{r)f% 2 \r) (12) 
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To control the size of the space, a number No is intro- 
duced and 2n + 2N + I + L < N Q is required for all 
the basis functions. There are two choices to obtain 
symmetrized eigenstates of Hi n . In the first choice, the 
set $j m rig(r, R) is firstly symmetrized and orthonormal- 
ized before carrying on the diagonalization. However, 
this procedure is complicated. Therefore we make the 
second choice, in which the set $j m nq(r, R) is simply 
used without symmetrization but with the requirement 
that I must be even and all the basis functions satisfying 
2n + 2N + 1 + L < N are included without missing. This 
requirement assures that the space is close under per- 
mutations, and no basis functions that will contribute 
to the symmetrized eigenstates would be missed, unless 
they are too high to have 2n + 2N + 1 + L > Nq. However, 
in this choice, a number of unphysical eigenstates with 
confused symmetry will emerge together with those with 
correct symmetry. Therefore a discrimination is needed 
as shown below. 

Let an eigenstate be denoted as ^ j m m where i is a 
serial number of the Jmll- series. Expanding in terms 
of the basis functions, 

Vj m m = £/™ n $j m n g (r, R) (13) 
i 

where the I in q must be even, and the coefficients B'/™ 11 
can be directly known from the diagonalization of Hin- 
If 4" j m Yii is correctly symmetrized, the coefficients would 
obey 

i 

(for all the q' with V even) (14) 

and 

Bt™ n A q J{p a) = 0, (for all the q' with I' odd) 

(15) 

In the summations of Eqs. (|T4")) and (ITSI) . I (in q) is re- 
stricted to be even. With the help of Eqs. (fT4")) and (H"5"1) . 
the states with confused symmetry can be discriminated 
and dropped, and all the symmetrized eigenstates under 
the restriction caused by Nq can be extracted without 
missing. They are one-to-one identical to those obtained 
via the first choice if the same No are used. In what 
follows \P j m m denotes only the symmetrized eigenstate, 
and the associated energy is denoted by Ejm. 



IV. EVOLUTION AND THE REPEATING 
3-BODY COLLISIONS 

With the eigenstates it is straight forward to obtain 
the time-dependent solution of H evo i as 

= e - iH °"° ' T */ = * c (R c> t)* <B (r,R J t) (16) 
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= Yl V J q r n (t) $ JmW (r,R) (18) 

where 

V J / lY1 (t) = Gj m .n q e- lTE '"' B^Bif 1 (19) 

i,q 

and r = uit. Obviously, Eqs. <p~6]> to (|T8|) give only an 
approximate solution because the set 4 j m m obtained via 
diagonalization each would deviate more or less from the 
corresponding exact eigenstate, and because only finite 
number of 4 j m m are used in the expansion. However, 
it is believed that, when the number of basis functions 
becomes larger and larger, the above ty(t) would be closer 
and closer to the exact solution. The crucial point is the 
convergency. In this paper the interaction is assumed to 
be weak. It turns out that, in this case, the convergency 
is satisfying as shown below. 

We shall demonstrate that the evolution is a repeating 
3-body collisions, and the effect of interaction will be 
also shown. For these purposes, we extract the following 
quantities from vE'(t). 

(i) The density p r (r,t) = J dfl r ^*(t)^(t), where the 
integration covers all the degrees of freedom except dr. 
Therefore J dr p r (r,t) = 1. Obviously, p r (r,t) is the 
probability density that the inter-distance is r. 

(ii) The density p R (R,t) fulfilling / dR p R (R,t) = 1 

(iii) The density p 9 {9, t) fulfilling / sin 9-d6 p e (6,t) = 1, 
where 6 is the angle between r and R. 

(iv) The density ((R,t) fulfilling J dR C(R,*) = 1. 

It was found that the cm. is distributed very close 
to the origin, therefore R is approximately proportional 
to r3. Therefore, the behavior of the particle 3 can be 
roughly understood via £(R, t). Incidentally, the behav- 
iors of all the particles are the same due to the sym- 
metrization. The analytical expression of these densities 
are given in the appendix. 

In order to have numerical results, as an example, the 
interaction is assumed to be a repulsive core as V{r) = Vq 
if r < 0.4, or zero if r > 0.4, where Vq is a constant to 
be given. The other parameters are chosen as r/ = 1.5, 
|a| = 1.5, and No = 20. To show the initial localization of 
the particles, we define the one-body density of the initial 
state as /5i(ri) = J dr2dr3 , &* 1 '$> j . This density is plotted 
in Fig.l. Starting from 4/, the details of evolution are 
given as follows. 

pu(R,t), p r (r,t), and pg(9,t) against t in the earliest 
stage of evolution is plotted in Fig2a to 2c. Where Vq = 
0.08 is assumed, and wt is from to 2ir (say, if u; = 
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FIG. 1: Contour diagram of the one-body density pi(ri) of 
the initial state plotted in the X-Z plane, rj = 1.5 and |a| = 
1.5 are adopted. Each maximum is marked with a x. The 
three particles form a regular triangle with side-length ~ 2.6 
initially. (Color online) 




FIG. 2: Contour plots of pn(R,t), p r (r,i), and pe(9,t) given 
in (a), (b), and (c), respectively, with rj — 1.5, |a| = 1.5, 
Vo = 0.08 and No = 20. Each maximum is marked with a x . 
(Color online) 

1000 x 2tt, then t is from to 0.001 sec). When t = 0, 
the inter-distances among the particles are about 2.6 as 
shown in Fig.l and 2b. When the evolution begins, the 
peaks of p r and pn move inward synchronously, while pg 
remains to peak at 90°. It implies a contraction of the 
regular triangle. Accordingly, the three particles rush 
towards the center leading to a 3-body head-on collision. 
When uit is close to ir/2, the inter-distances become much 
shorter, and pg becomes nearly uniform. It implies that 
the geometric character (i.e., the regular triangle) will be 
spoiled when the particles are close to each other. 

After the first collision, the particles move outward. 
When uit — tt, each particle will be close to the op- 
posite end of its initial position as shown by the solid 
curve of Fig. 3a. Where C(R., t) is given at the point 
R = (R,9r,(I)r) = (2.25, 7r, 0) and uit is given in the 
interval [0, 2tt\. The sharp peak of the solid curve demon- 
strates clearly that a particle arrives at the opposite end 
when tot ~ tt. Due to the symmetrization of the wave 



FIG. 3: C(R-i t) against t given in three selected intervals of 
t. R is fixed at (2.25, n,0) (i.e., R is lying along the negative 
Z-axis with R — 2.25. Meanwhile particle 3 is close to the op- 
posite end of the up- vertex of the initial triangle) . The three 
curves are associated with the three choices of A, namely, 0, 
200, and 400. When A = 0, the interval of cot is [0, 2tt] (solid 
curve). When A = 200, it is [200tt, 202tt] (dash curve). When 
A = 400, it is [400tt, 402tt] (dash-dot-dot curve). Vo = 0.08 
(a) and —0.08 (b). The other parameters are the same as in 
Fig. 2. (Color online) 



function, this is also true for other two particles. Then, 
the above process begins to reverse. When uit = 2ix, the 
system recovers its initial status nearly. If the interaction 
is removed, the recovery is complete and the system will 
undergo an exact periodic motion with the period 2-k/lo. 
The motion is characterized by the repeating head-on 
3-body collisions, each occurs once within the interval 
tt/lo. 

However, due to the interaction, the recovery is not 
exact. When the time goes on the effect of the weak 
interaction will accumulate and gradually emerge. This 
is shown by the dash-curve of Fig. 3a which describes the 
behavior of C(Rj t) after 200 rounds of head-on collisions. 
Where the peak is lower than that of the solid curve. It 
implies that the density is diffusing. Furthermore, the 
peak of the dash-curve has shifted a little right implying 
that the arrival is a little delayed. The diffusion and the 
delay will become more explicit when t is larger as shown 
by the dash-dot-dot curve. When the interaction is at- 
tractive, the above repeating 3-body collisions remain, 
and the densities remain to be diffusing. However, in- 
stead of a delay, the peak will arrive at the end earlier as 
shown in Fig. 3b, where Vq = —0.08. Due to the diffusion 
of the density as shown in Fig. 3, the phenomenon of re- 
peating 3-body collisions will become ambiguous when t 
is sufficiently large. 

When Vo is given at a number of values and R is fixed 
at (2.25, 7r, 0), C(R-7*) against uit is plotted in Fig. 4, in 
which uit varies in an interval close to 20l7r. In this figure 
both the locations and the heights of the peaks depend 
on the strength Vo. In general, a larger \Vq\ will cause 
a stronger diffusion and therefore a shorter life of the 
phenomenon of clear repeating 3-body collisions, and a 
more positive (negative) Vo will cause a larger shift of the 
peak to the right (left). In addition to Vo, the evolution 
depends also on the range of interaction. Therefore, by 
observing the time-dependent densities, information on 
the interaction can be obtained. 
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FIG. 4: t) against t in an interval close to 2Q1-k/uj. R 

is fixed at (2.25, 71", 0). Vo is given at four values marked by 
the curves. The other parameters are the same as in Fig.2. 
(Color online) 

The accuracy of the above numerical results depends 
on TVo, the number controlling the total number of basis 
functions. As an example the values of £(R, t) are listed 
in Table 1 to show the dependence. 



A crucial point in the device is the initial localization 
of the particles. If they are better localized (by increas- 
ing ujp) and/or they are more separated from each other 
initially (by increasing |a,-|), the repeating 3-body col- 
lisions would become more explicit and be maintained 
longer. Furthermore, if H evo i has a larger w, the inter- 
val between two successive collisions (~ 7r/w) becomes 
shorter, therefore the evolution would proceed swifter. 

Although the three particles are assumed to be iden- 
tical bosons, a generalization to fcrmions and/or distin- 
guishable particles is straight forward. In principle, the 
above device could also be used to study the three-body 
forces if they exist. 
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TABLE I: The values of C(R, t) with R = (2.25, tt, 0) and 
with three choices of iVo- The other parameters are the same 
as those for Fig.2 



No 


16 


20 


24 


C(R, 3l7r/w) 


0.0590 


0.0586 


0.0586 


C(R,9l7r/cj) 


0.0547 


0.0542 


0.0543 


C(R, 151tt/w) 


0.0472 


0.0468 


0.0469 



The convergency shown in this table is satisfying. Thus 
we conclude that the numerical results obtained by using 
Nq = 20 are accurate enough in qualitative sense. 

V. SUMMARY 

A model of a device containing repeating 3-body colli- 
sions is proposed to study the weak interactions among 
neutral atoms and molecules. The advantage is twofold, 
(i) The initial status can be precisely controlled, (ii) 
The weak effect of the interaction can be accumulated 
and therefore easier to be detected. Numerical results 
support that information on interaction can be thereby 
extracted. 



Appendix A: Analytical expressions of the 
time-dependent densities 

The analytical expressions of the time-dependent den- 
sities are as follows. 

(i) p R (R, t) = ^^(cNc) 2 ^ S n / n 6in5L'L 

jv c J,m,n,q',q 

x[v J q r u {t)]* v J q mU {t) 

xR*fW\R)fW L s) (R) (Al) 

(ii) Pr (r,t) = ^^2{cn c ) 2 ^2 5n>n&i>i5l>l 

N c J,m,Yl,q' ,q 

x[vp nI1 (t)]* v J q mYl {t) 
*r*ei 2 \r)fT\r) (A2) 

(iii) For the derivation of pe(9,t), a transformation is 
made so that r and R are transformed to r, R, 9, and 
the three Euler angles specifying the orientation of the 
triangle formed by the three particles. Then, integrating 
all the degrees of freedom except 9, we have 
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>*) = ^E( c ^) 2 E P^ mn '(*)r v j q mn {t) (-iy +l '- J 

N c J,m.n' ,II,q' ,q 

x[(2L' + l)(2L + l)(2r + l)(2Z + l)] 1 /2 [J r 2 dr f^, 2 \r) f% 2 \r) 
x[jR"dRf^}(R)f^ L 3 \R)} 



X 

A 



£ W(W'LL'; AJ) J^^C&voCftvMO, 0) (A3) 



where the Wigner and Clebsch-Gordan coefficients are 
introduced. 



J 



(iv)c(R,*) = ^E( c ^) 2 E 8 n , n sM^y u '(t)rv j q ^(t) 

N c J' ,J,m' ,m,W ,11, q' ,q 

^[Q^fc™ 'm'-kPlkX m-k^l' ,m' -k{R)YL,m-k{R)] 



X 

k 



xf%/?)(R) f^iR) (A4) 
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